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Abstract: We show that a couplet, a pair of closely spaced photon lines, in the X-ray spectrum is 
a distinctive feature of lepton flavored dark matter models for which the mass spectrum is dictated 
by Minimal Flavor Violation. In such a scenario, mass splittings between different dark matter flavors 
are determined by Standard Model Yukawa couplings and can naturally be small, allowing all three 
flavors to be long-lived and contribute to the observed abundance. Then, in the presence of a tiny 
source of flavor violation, heavier dark matter flavors can decay via a dipole transition on cosmological 
timescales, giving rise to three photon lines. The ratios of the line energies are completely determined 
in terms of the charged lepton masses, and constitute a firm prediction of this framework. For dark 
matter masses of order the weak scale, the couplet lies in the keV-MeV region, with a much weaker 
line in the eV-keV region. This scenario constitutes a potential explanation for the recent claim of 
the observation of a 3.5 keV line. The next generation of X-ray telescopes may have the necessary 
resolution to resolve the double line structure of such a couplet. 
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Haiku for a couplet 
better than one line 
is to find three, five, or nine 
we’d settle for two 


1 Introduction 

Increasingly precise cosmological measurements indicate that about 80% of the matter density of 
the universe is composed of constituents that are non-baryonic, and neutral under both color and 
electromagnetic interactions. However, the precise nature of the particles that make up this dark 
matter (DM) remains a mystery. One theoretically appealing possibility is that DM is composed of 
Weakly Interacting Massive Particles (WIMPs), particles with mass of order the weak scale that have 
interactions of weak scale strength with the standard model (SM) fields. This scenario is compelling 
because, provided the WIMPs were in thermal equilibrium with the SM at early times, just enough of 
them survive today as thermal relics to account for the observed dark matter density. 

While the WIMP framework requires that DM have interactions of weak scale strength with the 
SM fields, efforts to produce it at high energy colliders have proven fruitless so far. Likewise, efforts 
to directly detect DM in the laboratory through its scattering off nucleons, in spite of the increased 
sensitivity of current experiments, have also been unsuccessful. There are some tentative hints from 
indirect detection of DM annihilation to SM today, but there is no conclusive signal. In the wake of 
these null results, DM scenarios that retain the cosmological success of the WIMP framework while 
satisfying the current experimental bounds have become increasingly compelling. 

The matter fields of the SM (Q, U‘’, D”’, L and E^) are known to come in three copies, or flavors. 
Different flavors carry the same charges under the SM gauge groups, but have different couplings to 
the Higgs, and so differ in their masses. One interesting possibility, which has been receiving increased 
attention, is that DM, like the SM matter fields, also comes in three flavors [1-11] or has flavor- 
specific couplings to the SM [12-15]. Specific DM candidates of this type include sneutrino DM in 
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supersymmetric extensions of the SM, and Kaluza-Klein neutrino DM in theories with a universal 
extra dimension. 

In [2] , the simplest theories of flavored dark matter (FDM) were classified, and labelled as lepton 
flavored, quark flavored or internal flavored, based on the form of the interactions of the dark matter 
candidate with the SM fields. Models in which the DM has tree level interactions with the SM leptons 
but not with the quarks - as in lepton FDM - can naturally account for the observed abundance 
of DM while remaining consistent with all experimental bounds [16-18]. The reason is that strong 
production at a hadron collider or scattering off a nucleus rely on DM-quark interactions, which are 
loop suppressed in this scenario. In addition, because the average number of photons generated in 
DM annihilation to hadrons is much larger than in the case of DM annihilation to leptons, the limits 
from indirect DM searches using gamma rays are also weaker. Some other indirect signals of DM 
annihilation, such as the positron flux, are enhanced for lepton FDM, but regions of parameter space 
for which the DM is a thermal relic remain viable. 

In general the couplings of lepton FDM violate the flavor symmetries of the SM. In order to avoid 
conflict with the very stringent bounds on flavor violating processes such as ^ > ey, while giving 

rise to an annihilation cross section of weak scale strength, the couplings of DM to the SM leptons 
must be aligned with the SM Yukawa interactions. In theories where the flavor structure is consistent 
with Minimal Flavor Violation (MFV), so that the only sources of flavor violation are associated with 
the SM Yukawa couplings, this condition is automatically satished. Then each dark matter flavor is 
associated with a corresponding lepton flavor. It is this class of theories that we shall focus on. 

In realizations of FDM that respect MFV, the mass splitting between a pair of different DM flavors 
is dictated by the corresponding SM Yukawa couplings. In simple models, this splitting is proportional 
to the difference in the squares of the Yukawa couplings, so that for a Dirac fermion we obtain, 

^x,i - ^x,3 - V iVi - 2 /|) • ( 1 - 1 ) 

In this expression the index * = 1, 2, 3 runs over the three lepton flavors e, /r, r. While represents 
the mass of the DM flavor, yi represents the Yukawa coupling corresponding to the lepton 
flavor. The constant rj has the dimensions of mass and depends on the dynamics which UV completes 
flavor at some high scale A. If threshold effects at this scale generate mass splittings at tree level, rj 
can naturally be of order m^, where is the average DM mass. In this case the largest splitting, 
that between the e and r flavors, is expected to be in the MeV-GeV range for weak scale DM. If 
tree-level contributions at the threshold are absent, the leading effects are then loop suppressed. The 
largest splitting is then much smaller, in the keV-MeV range. 

Since the Yukawa couplings are small, the splittings between the different DM flavors are sup¬ 
pressed relative to the mass of each flavor. If the splittings are smaller than the electron mass, the 
dominant flavor-conserving decay mode 

Xi Xj + Vj (1.2) 

is slow on cosmological timescales, so that the lifetimes of the heavier flavors are much longer than the 
age of the universe. Then all three flavors are expected to contribute to the observed DM abundance. 

Now, suppose a tiny additional source of explicit flavor breaking is present in the theory, so that 
the flavor violating decays 

X^^Xj+l (1-3) 

can occur on cosmological timescales and dominate over the flavor-conserving decay. The monochro¬ 
matic photons produced in such decays then constitute a striking signal of DM. Provided this new 
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source of flavor violation is too small to contribute significantly to the splittings between the different 
DM flavors, the frequencies of the resulting gamma ray lines depend on the SM Yukawa couplings as 
in eq. (1.1). The constant of proportionality in eq. (1.1) cancels out when ratios of frequencies are 
considered. For example, if the r flavor of DM is the heaviest and Xe the lightest, we have, 

^ (Xr ^ Xu) ^ rnl-ml ^ 
w (Xr ^ Xe) rnl - ml 

Similarly, 

^(XM^Xe) rn-l-ml _ 

w(Xr-tXe) ml - ml ^ ml 

We see that this scenario predicts a pair of very closely spaced lines in the keV-MeV region corre¬ 
sponding to the Xt Xe and Xt x^^ transitions (the “couplet”), as well as an isolated line in the 
eV-keV region. Remarkably, the ratios of the frequencies of these lines are a firm prediction of this 
scenario. 

In the next section we review the MFV framework for models of lepton flavored DM. In section 3 
we choose a simple benchmark to illustrate the phenomenology of this scenario, focusing on constraints 
from various dark matter experiments and potential signals. We conclude in section 4. 

2 The Framework 

In this section we study the restrictions that MFV places on the parameters of theories of lepton FDM. 
The SM has an approximate C/(3)® flavor symmetry that acts on the three generations of fermions 
Q, C/'^, L and E^. This symmetry is explicitly broken by the SM Yukawa couplings. In extensions 
of the SM that respect MFV, any new parameter that violates the SM flavor symmetries must be 
aligned with the SM Yukawa couplings. Specifically, the Yukawa couplings of the SM are promoted to 
spurions that transform under the U{3)^ flavor symmetry. Any new interactions are then required to 
be invariant under this spurious symmetry. 

In the lepton sector of the SM there is an approximate U{Z)l x U(3)e flavor symmetry that acts 
on the left-handed SU{2) doublets and SU{2) singlets Ef. We denote U{2>)l indices by capital 
Latin letters and U{3)e indices by lowercase Latin letters. This symmetry is violated by the SM 
Yukawa interactions, 

Cy = Y^^L^H^E^+ h.c. (2.1) 

We can, however, make the Yukawa interactions formally invariant under the flavor symmetry by 
promoting the matrix * to be a spurion that transforms as (3, 3) under SU (3 )l x SU{3)e subgroup 
of U{3)l X U{3)e, and has appropriate charges under the U{1) factors. 

Theories of flavored DM posit a U(3)^ flavor symmetry that acts on the DM field x“. We use 
lowercase Greek indices to denote the DM flavor. We focus on the case where the DM field is a fermion 
that transforms as a singlet under the SM gauge interactions, and has renormalizable couplings to the 
to the SU{2) singlet leptons Ef through a scalar mediator (j). The mediator does not transform under 
the SM and DM flavor groups. The relevant interaction takes the form 

Cx = XJx^Efcl, + h.c. . ( 2 . 2 ) 

This interaction explicitly violates the U{3)e x U{3)^ symmetry. In general, it will give rise to lepton 
flavor violating processes at one loop. However, when MFV is imposed, the form of this interaction is 
restricted, with the result that flavor violating processes are forbidden. 
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Figure 1: a) Flavor preserving decay Xi b) A Feynman diagram contributing to the flavor 

violating decay Xi Xjl 


We impose MFV by identifying the DM flavor symmetry U(i)-^ with the U{3)e flavor symmetry 
that acts on the SU{2) singlet leptons in the SM*, and requiring that the new interaction be invariant 
under the spurious {7(3)® symmetry. This leads to a restriction on the form of = A/. To leading 
order in * we then have, 

A/=A(5/+AF;^y4' (2-3) 

MFV also governs the form of the DM mass matrix. If % is a Dirac fermion we can write the mass 
term in the Lagrangian as 

Cm = mJxpX'^ = m/xjX" ■ (2-4) 

MFV restricts the mass matrix to be of the form 

m/ = fhS/ + ^ (2.5) 

Then, in a basis where the lepton mass matrix is diagonal, we see that the splittings between the 
different DM mass eigenstates are governed by the SM Yukawa couplings in accordance with eq. (1.1). 


3 A Benchmark Model 


In this section we construct a simplified benchmark model which exhibits the features described above 
and consider its phenomenology in detail. We choose the dark matter y to be a Dirac fermion which 
transforms as a fundamental under the SM flavor group SU{3)e- The only additional terms in the 
Lagrangian beyond the SM are given by eqs. (2.2) and (2.4). We assume that the structure of the 
coupling and mass terms is restricted by MFV, and follows eqs. (2.3) and (2.5). For this benchmark 
we will further restrict to the special case where the DM mass terms and couplings generated at the 
high flavor scale A are universal, so that fh and A are zero at this scale. The mass and interaction 
terms for the DM in four-component notation at scale A are then given by. 


C = m^X^x'‘ 


^X*(l + 7®)eV'' +h.c. 


(3.1) 


Here e* is the four-component spinor corresponding to the charged leptons of the SM. The only free 
parameters in this model are the masses of the dark matter and the mediator, and a coupling A. As we 
shall see, the mass splittings generated in this case are finite at one loop and do not contain divergent 
pieces. Later we will introduce a tiny source of flavor violation. 

*One could also identify U{3);^ with the U{3)l symmetry that acts on the SU{2) doublet leptons (see [2]), but the 
main results do not depend on this choice. 
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Figure 2: One-loop correction to the two-point function of x 


In what follows we determine the splittings between the different dark matter flavors. We then 
compute expressions for the lifetimes of the heavier flavors in the presence of a small source of flavor 
violation. This model is constrained by a number of experiments. Constraints from g — 2 oi the 
muon and monophoton searches tend to be subdominant to direct detection constraints [18, 19]. LHC 
constraints on the production of the mediator (j) and indirect detection constraints from dark matter 
annihilations into positrons and photons can also be relevant. We study these signals in turn. 


3.1 Splitting and decays 

The breaking of the lepton flavor symmetry C/(3) —>■ C/(l)^ by the SM lepton Yukawa couplings is 
communicated to the DM sector through the FDM interaction, eq. (2.2). Even though they are 
assumed to be degenerate at tree level, mass splittings between the three x flavors will be induced at 
the loop level. In particular, let us consider the 2-point function for the three flavors of x (see figure 
2). The one-loop contributions are identical in the limit of massless leptons, but at order 0{m‘f) they 
begin to differ, thereby giving rise to differences in the wavefunction renormalizations for the Xi- Due 
to the chiral nature of the FDM coupling, it is easy to see that there is no direct mass renormalization. 
Once the helds Xi are brought back to canonical normalization, however, a mass splitting is induced 
between them. 


TOy i — 


. rriyX^ 

=Amu= ^ 


327r2 


dxx log 


X -I- (1 — x)m'j j — x(l — x)r 


xm'^ + (1 ~ x)m'j^ — x{l — x)r 


To leading order in the Yukawa couplings this yields. 


Amij 


647r2 


F{m\lm\), 


(3.2) 


(3.3) 


where yn denote the Yukawa couplings of each lepton and v = 246 GeV is the Higgs vacuum expectation 
value. The loop function F[x) is given by. 




(3.4) 


2 3 

We see that Xt is split significantly further from Xe and Xp. than these two states are split from each 
other. For an overall mass scale in the GeV regime and m^ ~ 0(100) GeV, 




(3.5) 

(3.6) 
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It is interesting to note that the one-loop splitting calculated above is a finite effect, suppressed by 
for large (j) masses. Note that the sign of Am is not arbitrary, and as a consequence Xt is the 
heaviest DM flavor. At two loops there arises a logarithmically divergent contribution to the mass 
splitting, where the log divergence is cut off at the UV flavor scale, A. This two-loop effect can become 
important for very large and A. We estimate that it is subleading to the finite one-loop calculation 
for the range of m^ we are interested in, provided the new physics scale A is less than 100 TeV. 

Once the x masses are split by these loop corrections, only the lightest x is truly stable. This is 
true even in the exact MFV limit where the D(l)^ flavor symmetry is preserved. For this benchmark, 
the splittings are smaller than the mass of the electron. Then, the leading contribution for flavor¬ 
preserving X decays arises at one loop and is illustrated in the left panel of figure 1. Note that this 
contribution is very suppressed due to the following three factors. 


• With a X mass splitting of order keV, the kinematically available phase space is extremely small. 
This results in a significant suppression for the 1 —> 3 process. 

• The loop amplitude is suppressed by the momentum-exchange scale, or more concretely by 
(Am/m^). 

• The lepton propagators in the loop couple to cj) on one end and to on the other. However, 
the former couples to right-chirality leptons while the latter couples to left-chirality leptons. 
Therefore both lepton propagators need a mass insertion to obtain a nonzero amplitude, so the 
decay rate is further suppressed by m^. ml.. 

As a consequence of these effects, heavier flavors are long-lived on cosmological time scales. 

Since the rates of flavor-preserving x decays are so extremely small, even a very small flavor- 
violating contribution can easily be the dominant channel for the decays of the heavier x- We add the 
following source of flavor violation to the Lagrangian, 

Cfv = ^SXij Xi(l -k -k h.c. (3.7) 


The leading mechanism for flavor-violating x decays are dipole transitions Xi Xi7) which are 
illustrated in the right panel of figure 1. Note that unlike the flavor-preserving decays, these are two- 
body decays, and there are no suppressions due to lepton masses. The rate for these flavor-violating 
decays can be calculated in a straightforward manner. If we assume that all off-diagonal couplings in 
dXij are of the same size dX <C A, then to leading order 


F-. = r 

^ ^ Xi- 


■Xil 


e^X^SX^ (Amy)^m^ 
1024 TT^ ml 


(3.8) 


where we have neglected higher order terms proportional to lepton masses. There are two aspects of 
this decay worth mentioning. First, 5X is a free parameter that allows the heavier flavors to decay, 
but is otherwise small enough to leave the phenomenology unchanged. Second, in order to not violate 
the parametric prediction of the line energies, should be small enough such that its contributions 
are subdominant to the previously calculated loop-induced splittings. This condition is not difficult 
to satisfy in specific cases. 

Note that F^ oc (Am^)^. Thus, the rate for the transition between x^ and Xe will be many orders 
of magnitude smaller than those between Xt and one of these two states (which have essentially the 
same rate as each other). Considering that the transition between x^i and Xe corresponds to both a 
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smaller energy and to a much smaller rate, we expect it will be practically unobservable. On the other 
hand, a robust prediction of this setup is that if an X-ray line signal is observed in the keV region, 
then at around 1% energy resolution (more precisely, an energy resolution of the order of yj^/y^) it 
will reveal itself as being made up of two line signals of comparable strength. 


3.2 Relic abundance 

The dark matter annihilation rate in each channel {xiXj ) is given by. 


(av) = 


X^mz 


327r(m2 -|- rri^) 


2\2 • 


(3.9) 


Since the dark matter is a Dirac fermion, there is no p-wave or chirality suppression, and thus the 
annihilation cross section today is the same as in the early universe to a good approximation. 

Since all flavor combinations of dark matter co-annihilate with one another with the same cross 
section, the cross section that gives rise to the correct relic abundance today is the same as for a single 
species of Dirac fermion DM, given by 


{av) = 2 X (2.2 X 10 cm^/s). 


(3.10) 


The factor of two relative to the canonical quoted value (for Majorana DM) arises due to the Dirac 
nature of the dark matter. The region of parameter space leading to the correct relic density is shown 
in figure 3 as a red band. 

Note that the above calculation applies only if we assume that the interaction with leptons alone 
is responsible for the dark matter thermal relic abundance. Any coupling to additional non-SM states 
will alter these numbers. This constraint can also be relaxed if the dark matter relic density is set by 
an asymmetry, which can arise rather naturally in these models [9]. 


3.3 Direct detection 


As discussed in detail in [2] , lepton flavored dark matter can scatter off nuclei via a one-loop photon 
exchange. These constraints can be severe in the region where the dark matter is a thermal relic. 
The dominant contribution to the WIMP-nucleon cross section is flavor diagonal and for each flavor 
of FDM it is given by. 






2 

1 + 3 log 



2 


(3.11) 


Here is the reduced mass of the dark matter-nucleon system and Kg represents the infrared cutoff 
in the loop calculation for the effective DM-photon coupling. This cutoff is mg, the mass of the 
corresponding lepton, unless mg is smaller than the momentum exchange in the process, of the order 
of 10 MeV. We therefore use K^ = to,- and A^ = m^, but set Ag = 10 MeV. In extracting constraints 
from the null results of direct detection experiments such as LUX, we use the total rate, summed over 
all three FDM flavors. The region of parameter space excluded by LUX is shown in figure 3 as the 
blue-shaded region. 


3.4 Indirect detection 

In the limit m-^ <C m^, (Am)^ and {av) both scale approximately as and therefore choosing 

a fixed mass splitting or requiring thermal relic abundance leads to potentially observable signals for 
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Figure 3: For mediator masses = 100, 200 and 400 GeV, we plot the position of the X-ray signal 
(in keV, gray dashed contours) as well as a number of constraints. Direct detection constraints from 
LUX are shown as the blue-shaded region, while the indirect detection constraints from positrons and 
photons are shown as the purple and green-shaded regions, respectively. The red band shows the 
region where correct relic abundance is obtained. 









indirect detection searches in photons and positrons (with the caveats mentioned at the end of the 
previous paragraph). The constraints from both indirect detection channels are more stringent for 
lower mass dark matter, since the signal rate scales as the square of the x number density, which 
itself scales as rn~^, while the background flux as a function of energy does not change as rapidly. 
Therefore, for a given Am or {(tv), these constraints can be weakened by increasing the dark matter 
mass and either decreasing the coupling or increasing m^. 

For the positron constraint from the AMS-02 experiment [20], the signal has contributions both 
from the prompt positrons produced when one of the annihilating particles is Xg, and also from 
secondary positrons from the decays of fj,'^ and that are produced when one of the annihilating 
particles is or Xt- The spectrum of the secondary positrons is shifted towards lower energies 
compared to the prompt positrons (not to mention that the branching ratio of r —>■ e + A is rather 
low), and therefore the bound from AMS-02 comes mostly from the prompt positrons. The bound is 
shown in figure 3 as the purple-shaded region. Note that the positron constraint is significantly weaker 
than the constraint from direct detection across the parameter space, and therefore the inclusion or 
non-inclusion of secondary positrons in determining the bound turns out to be academic. Owing to 
the relative factor of two between Dirac and Majorana DM (eq. (3.10)), the latter being relevant for 
SUSY for which the AMS bounds are calculated, the bound on the FDM annihilation cross section 
leading from prompt positrons (any one of the three xe flavors annihilating with Xe) is related to the 
total annihilation cross section (all nine annihilation channels) as 

{av) < 6{(Tv)bound,e+ , (3-12) 

where (cru)bound,e+ is the experimental bound quoted in [20] for a Majorana DM annihilating to e“'"e“ 
with 100% branching fraction. 

Similar to the case of positron constraints being most sensitive to prompt positrons in the final 
state, indirect detection in photons is most sensitive to r’s in the final state, since more photons 
are produced from r’s than from e’s or /r’s. One can therefore formulate the bound from indirect 
detection in the photon final state [21] in terms of the effective annihilation cross section leading to 
the production of r’s. 


(av) < 6{crv)bound,'f, (3.13) 

where {o'v)bound,-y is the experimental bound quote in [21] for a Majorana DM annihilating to t'^t~ 
with 100% branching ratio. The constraint from indirect detection in photons is shown in figure 3 as 
the green-shaded region. 

Dark matter annihilating to leptons can potentially have significant constraints from the CMB [22, 
23]. However, the annihilation into muons and taus has a low efficiency to inject energy into the 
CMB [24], and the constraints are subdominant to the other constraints considered above. 

3.5 The 3.5 keV line 

An X-ray line signal has been observed at 3.5 keV [25, 26]. We note that there is currently no 
consensus about the interpretation of this observation as arising from a dark matter signal [27-38]. 
Nevertheless, a large number of DM models have been proposed to explain this signal. The ideas 
that have been proposed include sterile neutrinos [39-64], axions [65-80], supersymmetry [81-86] and 
a number of other mechanisms [87-118]. We now show that the 3.5 keV line, if confirmed, can be 
easily accommodated within the framework of the minimal model. We show in figure 3 the region of 
parameter space that provides the right relic abundance for dark matter, consistent with a 3.5 keV 
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Figure 4: Constraints on the mass of dark matter, and the mediator, when the X-ray line 
is at 3.5 keV. The contours show the value of coupling A, and the red band shows the region where 
correct relic abundance is obtained. The blue and purple-shaded regions show the exclusion from LUX 
and from AMS respectively. 


line. Fixing the splitting to be 3.5 keV, we show the direct detection and the AMS positron constraint 
in figure 4 along with the region of parameter space consistent with the requirement of correct relic 
abundance. 

This scenario further predicts that closer inspection of such a line will reveal two closely spaced 
lines corresponding to the Xr Xfj, Xt Xe transitions. As noted in eq. (1.5), the ratio of the 
line energies in this couplet are set by the charged lepton masses. Therefore, 

Suj = Uj{Xr Xe) - ^{Xt Xfj) (3-14) 

TJl?' 

«a;o—= 12.4 eV (3.15) 

where wq = 3.5 keV is the average frequency of the two Xt decay lines. The broadening of the line 
due to the kinetic energy of the DM scales with its velocity. For typical astrophysical sources, the DM 
velocity ranges from (1-3)x 10“^, resulting in a broadening of 0(1-10) eV. Whether the double line 
feature gets washed out by the thermal broadening thus depends on the astrophysical source. While 
this splitting is not currently measurable, it is within the design resolution of upcoming experiments 
like ASTRO-H [119, 120]. The couplet constitutes a “smoking gun” signal of lepton FDM scenarios 
at these experiments. 

The lifetime for a decaying dark matter candidate to be consistent with the observed signal is 
given by (see for instance [93]), 

7 keV 

tdm — (10^® sec)-. (3.16) 

For = 150 GeV, one obtains tdm ~ 10^° s. For = 500 GeV and A ~ 1, this would require 
i5A ~ 10“®. The additional mass splittings introduced by this level of flavor violation are subdominant 
to the MFV contributions calculated above. 
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4 Conclusions 


We have studied models of lepton flavored dark matter within the MFV framework. In this scenario, 
mass splittings between different dark matter flavors can naturally be small enough that tree level 
decays of heavy flavors are kinematically forbidden. Then all three flavors of dark matter can be 
long-lived on cosmological time scales (the lightest flavor being exactly stable). The ratios of the mass 
splittings between the three possible pairings among the DM flavors are predicted. 

When even a very small source of flavor violation is present in the dark sector, a new decay channel 
becomes available for the decay of heavier dark matter flavors through a dipole transition. While the 
lifetime associated with this decay may be orders of magnitude longer than the age of the universe, 
it can still be the dominant decay channel and gives rise to a very distinct final state. These decays 
result in two very closely separated photon lines, the couplet. For weak scale DM, the overall energy 
of the couplet is naturally in the keV-MeV range, with the splitting of the two lines in the eV-keV 
range. 

We have focused on the detailed phenomenology of a specific model of lepton flavored DM. In this 
model, the mass splittings are radiatively generated by finite one loop effects arising from the breaking 
of the lepton flavor symmetry by Yukawa couplings. The sign of the contribution is fixed, with the 
result that Xr is the heaviest and Xe the lightest state. 

This scenario is a potential explanation for the claimed observation of a 3.5 keV line in the X-ray 
spectrum, and there exist regions in parameter space of the model where such an explanation is entirely 
consistent with the observed DM relic density as well as with experimental constraints set by a number 
of direct and indirect detection experiments. While a dark matter explanation of this line is in dispute, 
the smoking-gun double-line structure predicted by the couplet can be directly tested experimentally. 
In particular, the next generation experiments might be able to resolve any such feature in the X-ray 
spectrum. 
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